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Abstract. The generation of heralded pure Fock states via spontaneous parametric 
down conversion (PDC) relics on perfect photon-number correlations in the output 
modes. Correlations in any other degree of freedom, however, degrade the purity of 
the heralded state. In this paper, we investigate spectral entanglement between the two 
output modes of a periodically poled waveguide. With the intent of generating heralded 
1- and 2-photon Fock states, we expand the output state of the PDC to second order 
in photon number. We explore the effects of spectral filtering and inefficient detection, 
of the heralding mode, on the count rate, and purity of the heralded state, as 
well as the fidelity between the resulting state and an ideal Fock state. We find that 
filtering can decrease spectral correlations, however, at the expense of the count rate 
and increased photon-number mixedness in the heralded output state. As a physical 
example, we model a type II PP-KTP waveguide pumped by lasers at wavelengths 
of 400 nm, 788 nm and 1.93 (ira. The latter two allow the fulfillment of extended 
phase matching conditions in an attempt to eliminate spectral correlations in the PDC 
output state without the use of filtering, however, we find that even in these cases, 
some filtering is needed to achieve states of very high purity. 



PACS numbers: 42.65.Lm, 42.50.Dv, 42.65.Ky 
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1. Introduction 

Pure photon Fock states, in particular, single-photon states, are useful for many 
quantum-optical applications, including quantum information processing, quantum 
computing and quantum cryptography [TJ. Some novel uses for photon number states 
include the generation of other non-Gaussian states, such as Schrodinger kitten states 

A common method for creating single-photon states makes use of spontaneous 
parametric down conversion (SPDC), where a CW or pulsed laser is used to pump a 
nonlinear crystal. Due to the nonlinear properties of such a crystal, photons in the pump 
laser are spontaneously down converted into pairs of lower energy photons. Momentum 
conservation ensures that the photon wave packets are well localised with respect to each 
other. The detection of a single photon in one spatial mode (idler) heralds the presence of 
another single photon in the other spatial mode (signal). In practice, however, given high 
enough pump power, the presence of higher order photon-number terms in the output 
state can lead to a photon-number mixed state in the signal mode, when inefficient 
detectors mistake two (or more) photons for one. 

Energy conservation ensures that the frequencies of the downconverted photons 
always sum to the pump frequency. For CW pumped PDC, these correlations cannot 
be avoided, but pulsed pump light allows this constraint to be weakened. Strong 
spectral correlations are another potential source of mixedness - the signal state is 
projected into a spectrally mixed state when a frequency-insensitive detector heralds 
a single photon in the signal mode. In the context of single-mode versus multi-mode 
descriptions [U E] this property can also be interpreted as projecting the single photon 
onto different distinguishable broadband spectral modes [6]. In recent years, there 
has been a growing effort in engineering pulsed SPDC sources to produce photons 
uncorrelated in frequency, i.e. those with a separable joint spectral amplitude (JSA). 
Some examples include manipulating the crystal length, material, bandwidth and central 
frequency [3 El El EDI Ell El E21 E31 EH E5] as well as filtering the pump field, prior to 
down-conversion, using an optical cavity [IB] . Another promising technique produces a 
source of counter-propagating photons with a separable JSA [T7] . 

Typical theoretical analyses of multimode effects in PDC truncate the output state 
to first order in photon number. In this paper, we will extend our analysis to include 
second-order photon number contributions. We will compare the more humble method of 
pre-detection filtering, of which the main shortcoming is the loss in the production rate, 
with a more experimentally challenging approach where the JSA is made as separable 
as possible by satisfying various phase matching conditions. We find that even if these 
conditions are fulfilled, some level of filtering is still required in order to achieve states 
with very high purities. For both methods, we characterise the effects of higher order 
photon number contributions, on the generation of single-photon states. In addition, 
we consider the production of two-photon Fock states, conditional on the detection of 
two photons in the idler mode. 
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As a practical example, we consider a periodically poled KTP (PP-KTP) wave 
guide - this eliminates any spatial correlations - pumped by ultrafast optical pulses to 
ensure that the pairs exhibit tight localisations in time. 

This paper is structured as follows. In section |2j we introduce spectral notation for 
photon states as well as a theoretical model of spectral effects in SPDC. In section [3] we 
introduce spectral filtering. In section [4] we present analytical results for the probability 
of detecting a single photon in the heralding detector, the heralded output state, its g^ 
and purity, as well as the maximum fidelity between the heralded output state and an 
ideal pure state. In section [5j we present similar results for the generation of two-photon 
Fock states. In sections [6] and [7], we illustrate these results using realistic parameters. 
In section |8] we discuss our results. 

Finally, a note on nomenclature. In our theoretical calculations, we have a 
preference for using frequency (as opposed to wavelength) due to its direct relationship 
to energy conservation in SPDC, however, we have made an attempt to also present 
our results in nanometers (nm), for readers who are accustomed to "thinking in 
wavelengths" . All frequencies quoted in this paper are angular frequencies in units 
of s -1 . When describing Gaussian filters and beam profiles, we will specify the central 
frequency and Gaussian standard deviation in s _1 as well as the central wavelength and 
the FWHM in nm. 



2. Parametric Down Conversion 



A single-mode PDC can be modeled in the interaction picture, where the evolution of 
the state vector is given by [18] 

\ip(t)) = exp(Kt(a\al - aia s ))\ip(t )) . (1) 

This gives the output of a PDC [|J in the number basis, as 



j oo 

— — ^tanh(/rf)»» s , (2) 



COS— , , 

n=0 



where k is the effective nonlinearity and is a function of the pump power and the 
properties of the nonlinear crystal and t is the interaction time. The output state is 
correlated in photon number. 

To take the spectral properties of the system into consideration, we define the 
creation operator for a photon with a spectral distribution i/j(uj) as [5J [TU] : 

A\ = [ dw^{u)a\oj) . (3) 



A\p satisfies all the standard bosonic commutation relations, such as [A^jAt ] = 5k,k'i 
where and ipk'i^) are orthogonal spectral functions, i.e. J ipk(u)ipk'(u)dLj = <5 fcjfc /. 

| Note that a PDC state of one spatio-spectral mode is exactly equivalent to a two-mode squeezed 
beam. 
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An n-photon state can be written as 

l^> = 4t(^)1 > ■ ( 4 ) 



We emphasise the distinction between the states A^A^ = a/2|2;^) and 
A]p k Alp / = |1; ipk) |1; ipk') for k ^ k', where the former is a two-photon Fock state and 
the latter consists of two single-photon Fock states. 

For type-II down conversion, where the pump is non-delpleting, i.e. classical, we 
can take the multimode Hamiltonian to be 120] 



H{t) = [ d 3 r X {2) E(+\r,t)Et\r,t)Ei-\T,t)+R.c., (5) 
Jv 

where V is the spatial mode volume in the waveg uideand4(r,t) = i?j +) (r, t)+E^\r, t) 
are the three interacting fields with j — p,i,s denoting the pump, idler and signal modes 
respectively. 

E { p + \z, t) = Apf dcu p a(tu p )e i[k ^ )z+UJ " t] (6) 
M~\z,t) = I (ko j A{u j )at{u j )e- i ^^ )z+ "rt , (7) 



where a{ui p ) = exp(— {u p — fj Jp ) 2 /2a p ) is the pump envelope function and we have 
restricted the spatial integral to be over only one dimension, ie. z, and j = i, s. This 
Hamiltonian does not commute with itself at different times and therefore, the evolution 
of the state vector should be taken to be 

\m) = U(tM(t )) = Te-*tio dt ' H{t,) m o )) , (8) 

where T is the time-ordering operator. U(t) can be expanded into what is known as 
the Dyson series. This is very challenging, and as an approximation, we will drop the 
time-ordering operator and expand U(t) as a Taylor series. In the scenario presented 
in this paper, we do not expect this type of approximation to be problematic, however, 
we acknowledge that problems may arise when considering input states of a quantum 
nature, as has been investigated by Leung et al. plj. A description of multiple 
pair creation in degenerate SPDC has also been analysed in the Heisenberg picture 
by Wasilewki et. al [22] and Mauerer [23] • 

We will assume that A(ujj) is slowly varying over the frequencies of interest and 
therefore we can bring it outside of the integral. We can now write 

t poo p 

dt'Hit 1 ) = A dt' dz duidu s du p e- i[k ^ )+k ^ s) - kpM]z 

't J-oo J-L/2 J 

x e^+^-^a^aliu^iuJs) + H.c. , (9) 

where L is the length of the crystal and A = x^A p A(uJi)A(uj s ). For a pulsed laser, 
we can assume that the pump field, and therefore the interaction Hamiltonian, is zero 
before to and after t. Therefore we can extend the limits of the integration over time 
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to — oo and oo [24J. Performing the time integral yields 2n5(uji + uj s — uj p ) which then 
allows the u p integral to be evaluated, giving 

/ dt'H(t') = — 2nA / dz du>iduj s al(uji)al(uj s ) 

Jto J-L/2 J 

x a (wi + W8 ) e -<[M«<)+*.(«.)-M«i^.)]* + h.c. (10) 
Evaluating the integral over z yields 

jf dt'H(t') =2ixAL J duioluj s a\{ui)a\{u s )a{ui + u s )®{u h u s ) + H.c. (11) 

where $>(ui,u a ) is the phase-matching function 

/ LAk\ , , 

$(ui,u a ) = sine (— — , (12) 



V 2 

where sine (x) = sm(x)/x and Ak = ki(uJi)+k s (uj s )—kp(uJi+uj s ). For a periodically poled 
waveguide of periodicity A, Ak = ki{u)i) + k s {u) s ) —k p (uji+uj s ) + 2n/A [25]. Note that by 
picking the spatial integration to be centered around z — 0, it is possible to eliminate a 
phase term which would normally be present in equation (11). For the experiment, this 
corresponds to pre-chirping the pump pulse with an adapted phase progression. 

Following Grice and Walmsley [20], we Taylor expand the phase mismatch to first 
order such that Ak ~ Ak^ + k' s u s + k' i u i — k' p u p where Uj = luj — fj,j, kj = dkj{uj) / duj^^ 
and fij is the center frequency of a photon in mode j. We set /ij = /x s = \i and 
Hp = 2/i. We can achieve perfect phase-matching by picking A such that Ak^ = 
k s (n> s ) + ki(ni) — kp(fip) = 2n/A and therefore Ak w k' a v s + k[ui — k' p v p . 

To consider contributions from the 2-photon components of the down-converted 
state, we take the Taylor series expansion of the unitary evolution operator in equation 
(J8j) to second order (disregarding the time-ordering operator): 

U(t) « 1 + 1 / dt x H{h) + -^— 2 [ dt 2 H(t 2 ) [ dt 3 H(t 3 ) . (13) 

This gives the downconverted state 

I^pdc) =iv|(l + x 2 )|0)+x J J dudcuJ^^dl^dliu^O) 

+ Y J J duidu s f{ui,u s )d\{ui)d\{uj s ) 

(hj[du'J(ul,u' s )4(^M)\°)\ » ( 14 ) 




where x = 2nAL/ih and N is defined in equation (18). The joint spectral amplitude 
(JSA) is given by 

ffa, u s ) = N f a(u)i + u) s )$(u)i, u s ) , (15) 

where the normalisation parameter Nf is chosen such that J du>iduj s \f '(u^, uj s )\ 2 = 1. 

Any well-behaved complex function can always be decomposed in terms of a discrete 
basis of orthonormal functions (a well known example is the basis of Hermite functions). 
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This is known as the Schmidt decomposition. 

f(w 8 , u t ) = ^2 Mk(ui)( k (uj s ) , (16) 

k 

where the Schmidt modes and C&(u; s ) are normalised and may be complex and 

the Schmidt coefficients bk are real and \bk\ 2 = 1, if f(u a ,Ui) is normalised. It is 
useful to write the downconverted state in terms of the Schmidt decomposition (refer 
to Table [I] for creation operator definitions). 

|#dc) =n{(i + x 2 )|o) +xl>4 & 4jo) 

k 



.2 



k,k' 

where £fc(u^) are the Schmidt modes for the idler state and Cfc(^s) are the Schmidt modes 
for the signal state and 

-1/2 



N 



{|i + x 2 l 2 + lxl 2 + lxl 4 ( £ IW 2 + J>| 4 )} ■ (18) 



k,k' 
k<k' 



Notice, in the four-photon term of equation (17), when two photons are created in the 
same spectral mode (i.e. k = k') there will be a factor of v2 in front of each two-photon 
Fock state, increasing the probability of down conversion into such a state. This can be 



understood due to stimulation effects in the PDC process itself. Equation (17) can also 
be written as 



|#dc) =Ar{(l + X 2 )|0>+x$> fe l 1 ^l 1 ^ 

k 



+ E ^fc / |i;6)*|i;^fe / )il 1 ;Cfc)s|i;Cfc / )sj} ■ ( 19 ) 

k,k' 
k<k' 

We can characterise the spectral entanglement of the JSA by using the entropy of 
entanglement [26]. The entropy of entanglement can be defined, for the bi-partite state 

|*>=536 Jfc |l;&)|l;C*> (20) 

k 

in terms of the Schmidt values: 

Em) = -E & * lo s^) • ( 21 ) 

k 

The entropy of entanglement is valid only for pure bipartite states and, when defined 
in terms of the Schmidt decomposition, can not be applied to the entire output state in 



equation ( 19 ). However, we can apply it to the two-photon term to get some information 
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Creation Spectral State Description 
Operator Mode 



C T( k 


C*(^i) 


1; Ck)i 


initial idler Schmidt modes 




1; 6c)s 


initial signal Schmidt modes 


T(wi)( k (wi) 


|i;^C*>< 


filtered idler modes 


Dlt 

C 4>i 


R(uJi)( k (uji) 


1; R(k)i 


reflected filtered idler modes 


(f)j{Ui) 


|i; 


orthog. filtered idler modes 




<Pj(Ui) 




orthog. reflected filtered idler modes 




T~m{uJ s ) 


1 1 j 7~m ) s 


diag. single-photon signal modes 



Table 1. Summary of multi-mode creation operators, spectral modes and states. 

about the spectral entanglement arising only from the JSA. The entropy of entanglement 
ranges from zero for a product state to log 2 N for a maximally entangled state of two N- 
state particles, which in our case corresponds to a state containing N orthogonal spectral 
modes. In the limit of a maximally entangled JSA, the entropy of entanglement would 
be oo. 



3. Spectral Filtering 

A spectral filter can be modeled as a frequency dependent beam-splitter: 

d\uj) -> T(uj)c\uj) + R(uj)d\uj) , (22) 

where |T(a;)| 2 and |i?(c<j)| 2 are the transmitted and reflected probabilities and |T(c<j)| 2 + 
|i?(cj)| 2 = 1. In addition to the filter, we consider an inefficient detector which we model 
by a beam splitter of reflectivity 1 — r], followed by a perfect detector (refer to figure 
[T| . If the reflected mode of the filter and the reflected mode of the beamsplitter are to 
be traced out, the filter-beamsplitter combination can be modeled by a filter with the 
following transformation: 

d\oj) -> T(uj)c\uj) + R(uj)d\uj) , (23) 



where T(u) = T(Lu)y/rj and R(u>) — yl — \T{u)\ 2r q. In terms of the mode functions Cki 
this can be written as 

A^ T A + R A (24) 

where we have defined 



T Cfc = \ du}\T(u)( k (u)\ 2 (25) 



R Ck = J j dw\R(")(k(u)\ 2 (26) 

=-±- f dwT(u)( k (co)$(u) (27) 
Cfc 
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The definitions in equations (25)- (28) ensure that the creation operators for the filtered 
modes satisfy the commutation relations [Cr( k , CroJ = 1 an< ^ [-^flCk> -^J?cJ = ^ 
However, the filtered functions T( k (u s ) no longer define proper modes because the 
functions ( k (u s )T(u! s ) are, in general, not orthogonal to the functions ( k r(u s )T(u) s ) for 
k 7^ k! and therefore need to be orthogonalised (eg. using the Gram-Schmidt procedure) 
such that: 

3 
3 

where <pj{^i) are now the new modes defining the idler state and <fj{uJi) are the reflected 
modes that will be traced out, and 

u kj = J dw(f>j{w)*C, k {u)T{u) = T Cfc (l;0 3 -|l;TCfc) (31) 

v k j = J duipj(uY( k (uj)R(u) = R Ck (l; ipj\l; R( h ) . (32) 
The filter relationship in equation (plj) can now be written as follows: 



( 33 ) 

3 

Because the filter has been modeled as a frequency dependent beamsplitter, it will have 
similar properties to a beamsplitter. One property, that we will evoke throughout this 
paper, is the tendency to allow one photon, from an incident two-photon state, to pass 
through the filter while rejecting the other. Detection of a filtered two-photon state may 
then result in, even, a perfect detector mistaking it for a one-photon state. 



4. Generating single-photon Fock states 

Detection of a single photon in the idler mode heralds the presence of a single photon 
in the signal mode. In this section, we consider a hypothetical perfect detector, an 
inefficient detector and a Gaussian spectral filter placed in front of the lossy detector. 



4-1. Case 1: Perfect Detection in the triggering idler mode 

The POVM for a detector that perfectly distinguishes photon number, but gains no 
information about the frequency of the photon can be written as follows: 

n 1= I ^a t H|0)(0|aH = ^|l;0)(l;Ol- (34) 

J 3 

We can interpret this as: the detection of a single photon cJ(uj)\0), however due to the 
lack of spectral knowledge, u must be integrated over; or alternatively, the detection 
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(a) (c)(i) 



a(u) p ) 



/( 2 ) 



a+(w s ) 



o|(wi) 



Pfiit 



Pineff 





nonlinear crystal — perfect detector beam splitter ^^^^^ 



filters 



Figure 1. Schematic diagrams of SPDC setup for: (a) perfect detection in the 



triggering mode (see sections 4.1 and 5.1); (b) inefficient detection in the triggering 



mode (see sections |4.2| and 5.2); and (c)(i) filtering the triggering mode prior to 



detection with inefficient detector and (ii) the equivalent setup where the filter and 



beamsplitter have been combined into one filter (see sections 4.3 and 5.3) 



of a single photon in the spectral mode however due to the lack of knowledge 

about which mode it was in, it is necessary to sum over j. Since photon detection is 
destructive, the detected mode must be traced out. The probability of detecting a single 
photon in the idler mode using a perfect single-photon detector (refer to figure fTVa)), is 

pi = (^PDclIIil^PDc) = |Af |x| 2 • (35) 
Given a single-photon detection in the idler mode, the heralded signal state is 



- 1 Br i [ni|*PDc)<*PDc|n 1 ] = ^|6 fc | 2 |l;efc) s (l;efc|- (36) 



Pi 

Pi 

The g^ of the signal state, which we define as 

gm = ^ i A k t, j> (37) 

(^■(44)) 

is g^ =0. This reveals that there is only one photon in the signal mode, but not how 
pure it is. The purity of the heralded state is 

P 1 = Tr[pl} = \N 1 \ 4 Y,h\ 4 . (38) 

k 

For a state which only contains single photons, i.e. one that is heralded by a perfect 
detector, the purity is equivalent to the Hong-Ou-Mandel visibility 
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Without loss of generality, we order the Schmidt coefficients in decreasing order 
from k = 0. Therefore, the pure single photon state with the highest overlap with the 
projected state will be the photon mode corresponding to the highest Schmidt coefficent, 
60, and hence the maximum fidelity with a single photon Fock state is 

F 1 =maxF(p 1 ,\l;£ j )) = (1; fo| Pl |l; f > = N 2 • (39) 
j 

The fidelity does not depend on \- Increasing the strength of the nonlinearity only has 
an effect on how often the detector registers a click, however, once that happens, the 
signal mode is always projected into the same state. 



4-2. Case 2: Inefficient Detection in the triggering idler mode 

An inefficient detector can be modeled by the transformation in equation (J23j) where 
T(uj) = y/fj and R(u) = y/1 — rj, followed by a perfect detector (refer to figure pjb)). 
After the beam splitter, the joint signal-idler state is 

Pincff = Tr^[|^ incff ) (^ineffl] , (40) 

where 

pincff) = n{(1 + X 2 )|0) +xJ2 b ^Vv\Cl + y^\Dl)\0) 



+ y(E h ^'\\ {vclfil + (i - v)D\ k Dl 

k,k' 



+ VvV^liClDl + ClDl)) |0>) } . (41) 

The probability of detecting a single photon in the idler mode, using an inefficient 
single-photon detector, is 

Pi,ineff = Tr[Ilip inefr ] (42) 

= |iV| 2 | X | 2 r / {l + 2| X | 2 (l-r ? )( \hb«\ 2 + ■ ( 43 ) 



k,k' 
k<k' 



Given a single-photon detection in the idler mode, the heralded signal state is 

Pl.incff = " Tr^piPiaefflli] (44) 

Pl.incfT 

= |^V"l,ineflf| 2 | |^fc| 2 |l; £k)i (1) £fc| 

k 

+ 2| X | 2 (l-r ? )(^|6 fc | 4 |2;a) J (2;a| 
k 

+ Yl l^'l 2 ! 1 ;^! 1 ;^'}^ 1 ;^!^ 1 ;^'!)}. ( 45 ) 



k<k' 



where A^i )ino fj = Nx\fn/ yJPi,me&- The g( 2 > for this state, defined in equation (37), is 
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4M 2 (i-^)( E iw 2 +En 4 )- ( 4? ) 



k,k' 
k<k' 



The purity of the signal state is 

A,mefr = Tr[p 2 incff ] (4f 
= |iV det , liineff | 4 {E|6 fc | 4 + 4| X | 4 (l-r / ) 2 



x 



( E iw 4 +En 8 )}- ( 49 ) 



fc<fc' 



The maximum fidelity, between the heralded state and a pure Fock state |1; is 

F liineff = maxF(|l;^),p liincff ) = (1; £ |Pi,ineff|l; fo) = |iVi iine ff| 2 |6o| 2 (50) 



j 



4-3. Case 3: Filtering the idler state 

We now introduce a filter in the idler mode as shown in figure [l|c) . Applying a filter, 



as defined in equation (24), to the signal mode, gives the filtered state 



Pm = Tri>[\*at)(*at\], (51) 

where 

|* filt ) = N{ (1 + x 2 ) |0> + x E b kA{ k (u kj Cl. + v kj DQ |0> 



k,k',j,j' 

+ u kl v k , r Clpl 3i + W>^k) 1°) } > ( 52 ) 



and u k j and v k j are defined as per equations (31) and (32). The probability of detecting 
a single photon in the filtered idler mode is 

Pifiit = Tr[riip fi i t ] (53) 
+ E \bkbk'\ 2 {T 2 k R 2 y + T^ k/ R^ k + T kk >R k > k + T fe / fc R fcfe / j j , (54) 



fc,fc' 
fc<fc' 



where 



Tfcfc' = T Cfc T^,(l;TCfc'|l;TCfc) = E M ^'i ( 55 ) 

j 

R fcV = R <k R* Ckl {l;R( k ,\l;R( k ) = Y, v ti v h ■ ( 56 ) 
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Given a single-photon detection in the idler mode, the heralded signal state is 

Pi,m = — ^— TrcfniPfiitlli] (57) 
Pi,mt 

= l^ l!fil ti 2 {^M|T^|l;6)(l;^l 

k,k 

+ \X\ 2 ( E 6 fe 6^6|6|,T^R fc4 ititj 0)(0 |i e _i 5 _^ (58) 

where iVi^it = ^V"x/-v/Pi,fiit- Note that filtering the idler mode also changes the mode 
structure of the heralded signal state. The g^ for this state, defined in equation (37), 
is 

G (2) = 71 (Kq) 

9 ~ \Ni,m\ 2 (E k \hm k \* + lf y> [ > 

where 

lf = 2\ X \ 2 {2j2\bk\ 4 RkkTkk 

k 

+ E |&fe^fc'| 2 (RfefcT fe / fc ' + Rk'kTkk' + Rfcfc'Tfc/fc + Rfc'fc'T fcfc )| • (60) 

k < fc' 

The purity of the heralded state is 

P ljfilt = Tr[p 2 >filt ] (61) 
= l^,mt| 4 {El^r|T fc ,| 2 + | X | 4 {E2|6,| 4 



E l&fc&fc'HllfcfcT^ + R fc ,fcT fc ^| 2 + ^2|6 fc | 4 |R^T^ 

k, k' k 
k <k' 

+ E i & i & i'i 2 (E 2 i^i 4 i R ^ T ^+ R ^ T ^i 2 

+ E i & * 6 fc'i 2 i R fcfe T fc'ife'+ R fe'fc T fcfc'+ R feit' T ^ ( 62 ) 



k.k 1 
k <k' 



The density matrix in equation (58) is not diagonal in the |l;£fc) basis, but this can 
be easily achieved, for the part of the state which is relevant for calculating the fidelity 
with a single-photon Fock state, giving 

Pi,mt, P art = |A^i imt | 2 E & fc^ T fefcl 1 ;6)(l;4l = ^2d m \l;T m )(l;T m \ , (63) 



k,k 

where 



|l;r m ) = E^l 1 ;^)- (64) 
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and T m ((jj s ) are the new orthogonal modes defining the signal state. The maximum 
fidelity, between the heralded state and a pure single-photon state |1; 77), is 

Fi jfllt = maxFfll; n),pi >m ) = max (1; r z |pi. filt)Part |l; n) = maxd m . (65) 

l l in 

In an experiment, the spectral distribution r m should be chosen in any interferometric 
experiment to optimise for the best performance of the heralded single photons. 

In the extreme case where T(u) = 5(u—fi), i.e. the filter picks out a single frequency 
/i, the fidelity tends to unity and the signal state tends to the pure state 

|*i,filt,<5) = A^fiit^&fcCfcG")! 1 ;^) = Vd~^\l;T m ) , (66) 

k 

as x ~ y 0. This implies that it is possible to obtain arbitrarily pure single photon states, 
with the use of spectral filtering and by ensuring the nonlinearity strength is low. 



5. Generating 2-photon Fock states 

In addition to creating single-photon states, it is becoming increasingly desirable to 
create higher photon-number Fock states. In this section we will investigate the effects 
of detector efficiency, and filtering of the idler mode, on the generation of two-photon 
Fock states in the signal mode conditional on the detection of heralded two-photon 
states in the idler mode. 



5.1. Case 1: Perfect Detection in the triggering idler mode 

The projector for detecting two photons in any spectral modes (j and Qi will be 
separated into two parts: the part which detects two photons in orthogonal modes 
and the part which detects two photons in the same mode: 

n 2 = X>;G)(2;0I+ E |i;0)|i;0)(i;0l(i;0l- (67) 

3<f 

Refer to the schematic in figure [IJa). The probability of detecting two photons in the 
idler mode, with a frequency insensitive detector, is 

P2 = (*PDc|n 2 |vi>p D c) = |Af |xl 4 ( IW 2 + E l 6fc l 4 ) • ( 68 ) 

k,k' k 
k<k' 

Given a two-photon detection in the idler mode, the heralded state in the signal mode 
is 

p 2 = -TrJiI^pDc) (^ PDC |n 2 ] (69) 

P2 

= |iV 2 | 2 {]T|& fc | 4 |2;£ fc > s (2;£ fc | 

k 

+ \hh>\ 2 \l;Z k )s\i;Zk>)s{l;Zk\{l;&\} (70) 

k<k' 
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\N, 



h\ 2 At |0) (0\A S(k ® \h'\ 2 A^ |0) (0|i , (71) 



12 

2~ ^ 

fe fc' 

where N 2 = Nx 2 / y/p~2- The purity is 

p 2 = Ti[ P 2 2 ] = \N 2 \ i { Md'+Ew 8 }- ( 72 ) 

fc,Jfe' fe 

fc<fc' 

The g^ 2 -* for this state is g^ = 1/2. It is interesting to note that the g^ does not 
depend on the purity of the two-photon state. It will always remain at the value of 1/2 
regardless of whether the two-photon state is in a Fock state or in some other form. 
This reflects the fact that is only sensitive to the photon number, but not the modal 
properties of the state. 

The maximum fidelity between the heralded state and an ideal two-photon Fock 
state 1 2; £j) is 

F 2 = maxF(p 2 ,|2;^» = (2; f |p 2 |2; &> = l^ 2 | 2 |6 | 4 - (73) 

3 

5.2. Case 2: Inefficient Detection in the triggering idler mode 

The probability of detecting two photons in the idler mode, with an inefficient frequency 
insensitive detector (refer to figure [jjb)), is 

p 2 , incff = Tr[n 2 p ineff ] = |Af | X |V ( I + Yl l 6fc ! 4 ) • ( 74 ) 

k,k' k 
k<k' 



where pi ne fr is defined in equation (44). Because our analysis only extends to second 



order in photon-number, the expressions for the fidelity and purity will be the same as 



they were in section 5.1, where a perfect detector was used. If we included higher order 
terms, we would expect the fidelity and purity to vary as a function of \ and rj in a 
similar fashion to the single-photon case in section |4} 

5.3. Case 3: Filtering of the idler state 

We now introduce a filter in the idler mode as shown in figure [Tj(c). After filtering the 
state, the probability of detecting two photons in the idler mode, with an inefficient 
frequency insensitive detector, is 

P2,fiit = Tr[n 2 pmt] (75) 

= |iV| 2 |xl 4 (E^I 4 H T cJ 4 

k> 

+ l^'| 2 (I^J 2 |T Cfc ,| 2 + T^T fc , fc )), (76) 

k,k' 
k<k' 



where p^t is defined in equation (51). Given a two-photon detection in the idler mode, 



the heralded state in the signal mode is 
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P2,m = — Trc[n 2 pmtn 2 ] (77) 

P2 



|iV 2 ,mt| 2 {^^fT^|2;a)(2;^ 



+ ^hh'bfT k ~ k T k r k \l; 6) 1 1; £*) (2; fc| 
fc<fc' 

+ ^^^,T^,|2;e fe )(l;4l(l;^l (7 

+ X] ^^(T fc ^ + T fc ,^)|i;e fc )|i;ev)(i;4l(i;^ 



k,k',k,k'_ 
k<k',k<k' 



»2,filt 



— E^! T a i ki°)(°l i ^ 



Ar.fc 



® X) ^' & I' T fc',fc'^L, |0> <0| A , (79) 



fc' k' 



where iV^mt = Ny 2 / ^Jp2,m- Again g^ 2 ' = 1/2 as there are always two photons in the 
state. The purity of the heralded state is 

PW = Trfe,t] = |AW{ E l fo ^'i 2 (E 2 l^| 4 | T a T ^| 2 



fc,fc' 



~t~ / l^^' I I Tr, juTt,/ ?,/ -(- Tr,/ ilTt 7/ I 

+ Ei^i 4 (Ei^i 4 i t ^i 4 + E 2 i^i 2 i T M T ^i 2 )}- ( 8 °) 



k,k/ 
k<k' 



The density matrix in equation (78) is not diagonal in the |l;£fc) basis. This can be 
easily achieved, giving 

P2,mt = W d m d m >Ai m Ai |0> (0\A STm A ST , (81) 

= ^<|2;r OT )(2;r m | (82) 

+ d m d m /|l;r m )|l;r m /)(l;r m |(l;r m /| , (83) 



m, m 
m<m' 



where |l;r m ) is defined in equation (64). The maximum fidelity, between the heralded 
state and a pure two-photon state with an optimised spectral distribution function |2; tj) , 
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is 

^2,mt = maxF(|2; 7}),p 2 ,fiit) = max (2; Tj|p 2 ,fiit|2; n) = maxd 2 m . (84) 

( l m 

In the extreme case where T(u) = 5(uj — fi), i.e. the filter picks out a single frequency 
(jl, the fidelity tends to unity and the signal state tends to the pure state 

k k' 

= d m \2; T m ) , (86) 

as x ~^ 0- This implies that it is possible to obtain arbitrarily pure two-photon states, 
with the use of spectral filtering and by ensuring the nonlinear ity strength is low. 



6. Physical example I - correlated JSA 

As a physical example, we model a type II PP-KTP waveguide of length L = 3.6 mm 
and a periodicity of A = 8.8 /mi, pumped with a 400 nm laser with a lnm FWHM 
(op = 5.00xl0 12 s _1 ) which down converts to 800 nm in the signal and idler modes. 
In figure [2] (a)-(c), we have plotted the pump function, the PMF and the JSA for the 
given parameters. Figure [2] (d) shows the corresponding Schmidt modes: initial signal 
modes and initial idler modes. In addition, it can be seen that after filtering the idler 
state, the idler Schmidt modes take on different spectral shapes. These filtered idler 
modes are no longer orthogonal to each other and therefore need to be orthogonalised 
giving the orthogonalised idler modes. When the idler mode is detected, the signal 
state gets projected into a mixture of orthogonal modes, as shown by the diagonalised 
signal modes. Filtering and detection of the idler state changes the spectral shape of 
the signal state even though there is no physical interaction. This is a typical effect of 
entanglement. 

In general, the Schmidt decomposition can not be found analytically, but can be 
calculated numerically by computing the singular value decomposition of a discretised 
JSA. Unless stated otherwise, the results in this section were obtained using an 800 x 800 
grid, ranging over 0.2 x 10 15 s _1 , centered around Ui = uo s = //. We note that an 
insufficiently fine grid, or insufficiently large region, will result in inflated values for the 
purity. The entropy of entanglement for this particular JSA is E = 4.6. 

In this section we present results for: the probability of detecting a single photon 
in the idler mode; the g^ and purity of the heralded state in the signal mode; and 
the fidelity between the signal state and the desired ideal Fock state. We compare 
results for: an unfiltered idler state; an idler state filtered with a Gaussian filter 
T{uii) = exp(— (fif — djj)/2cr|), of various widths 07 and centered at the central idler 
frequency, where the filter function has been scaled such that the maximum value is 
always 1; as well as the limiting case where T(u>i) = 8(uJi — fif)- We also present similar 
results for heralding a two-photon state conditional on the detection of two photons in 
the idler mode. 
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Figure 2. (a) Gaussian pump function a(u>i + uj s ) with a lnm FWHM at \x v = 400 
nm. (b) Phase matching function <5>(wi,aj s ) for waveguide of length L — 3.6 mm and a 
periodicity of A = 8.8 /im. (c) The resulting JSA f(u>i,uj s ) = a(u>i + uj s )^(uji 1 uj s ). The 
JSA has been plotted as a function of the frequency, however corresponding values for 
the wavelength have been included, (d) Schmidt numbers and modes for the JSA (top 
to bottom): the first 20 Schmidt numbers bk', the first 5 Schmidt modes £fc(w s ) for the 
signal state; the first 5 Schmidt modes Cfc( w i) f° r the ilder state, as well as a Gaussian 
filter function of width ay = 2xl0 12 s _1 (dashed line); the filtered Schmidt modes 
T(uJi)(k(ui) for the idler state; the othogonalised idler modes (f>j(u)i); the diagonalised 
signal modes T m (u s ) 



6.0.1. Generating single-photon Fock states Due to the second-order truncation of 
the PDC output state, we are not considering 6- (or higher) photon contributions. 
At x — 0.5, the fraction of 6-photon states, to 2-photon states, would roughly be 
xVx 2 = 1/16. We will not plot results beyond that. 

In figure [3] (a), we have plotted the probability of detecting a single photon in the 
idler mode as a function of the nonlinearity \ an d the efficiency of the detector rj. Notice 
that the probability of detecting a single photon in the idler mode increases with higher 
detector efficiency and higher nonlinearity strength, as expected. 

The fidelity has been plotted in figure [3] (b). The inclusion of a filter has a drastic 
effect on the fidelity. It has a greater dependence on the strength of the nonlinearity, 
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than in the unfiltered case, however the overall fidelity is much higher. Notice as well 
that there is a trade off between the fidelity and the probability of detection. 

In figure [3] (c), we have plotted the g^ for the heralded state in the signal mode. 
Note that, for visual clarity, the figure orientation has been rotated by 7r around the 
z— axis, with respect to the other plots. It is useful to know that g( 2 > = for a single- 
photon state and 1/2 for a two-photon state. A curious thing is that decreasing the filter 
width results in higher fidelities, despite the higher proportion of two-photon states, as 
shown by the g^ 2 \ This suggests that for this particular JSA, the dominant cause of 
impurity is the spectral entanglement, rather than the resulting photon- number mixture 
due to the presence of higher-order terms. 

In figure |3] (d), we have plotted the purity of the state in the signal mode. Due to 
the 4-fold summation, the purity for the filtered case is very demanding computationally, 
therefore, we have only included examples of two filter widths. Results for the filtered 
case were computed using a 600 x 600 grid, ranging over 0.16 x 10 15 s _1 , centered around 
uji = u s = fi, and truncating bk with values below 10~ 2 . 

To achieve a fidelity of F = 0.95, using a heralding detector with efficiency 77 = 0.5, 
we could choose from a range of filter widths at different nonlinearity strengths. Different 



(a) (b) 




Figure 3. (a) The probability of detecting a single photon in the idler mode for (top 
to bottom): no filter; 07 = 3xl0 12 s _1 ; 07 = 2xl0 12 s _1 ; 07 = lxlO 12 s _1 . Note 
that the probability is plotted on a log scale, (b) The fidelity of the signal state with 
an ideal Fock state for (top to bottom): 07 = 0; 07 = lxlO 12 s _1 ; 07 = 2xl0 12 s _1 ; 
Of = 3xl0 12 s _1 ; no filter, (c) The <?' 2 ' of the signal state for (top to bottom): 07 = 0; 
07 = lxlO 12 s _1 ; 07 = 2xl0 12 s _1 ; 07 = 3xl0 12 s _1 ; no filter. Note the change in 
axis orientation, (d) The purity of the signal state for (top to bottom): 07 = 0; 
a j = lxlO 12 s" 1 ; 07 = 2xl0 12 s" 1 ; no filter. 
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combinations, however, result in slightly different probabilities of success. Figure [4] (a) 
shows the probability of success, and required nonlinearity, for a number of filter widths. 



(a) 
0.24 

0.2 

0.16 

* 0.12 

0.08 

0.04 



-2.4 

-2. 

- 1.6 f> 

S3 

-1.2 S 
0.8 
-0.4 



5 0.20 C 

a, (xl0 l2 .s--') 



(b) 
0.21 
0.2 
0.19 
0.18 
0.17 

0.16: 



0.5 0.6 0.7 0.8 0.9 
o- ( (*10 12 .r 1 ) 



-1.1 
- I. 

0.9 
-0.8 
-0.7 
Jo.6 



(c) 
0.23 
0.22 
0.21 - 
0.2 - 
0.19 
0.18 - 



0.4 



Figure 4. Probability of detecting a single photon in the idler mode, and required 
nonlinearity, for a number of filter widths, in order to achieve a fidelity of F — 0.95, 
using a heralding detector with efficiency r/ = 0.5, for: (a) correlated JSA (see section 
|6j; (b) symmetric JSA (see section 7.1|; (c) asymmetric JSA (see section 7.2|. 



6.0.2. Generating 2-photon Fock states At x — 0.25, the fraction of 6-photon states, 
to 4-photon states, would roughly be x 6 /x 4 = 1/16. We will not plot results beyond 
that. Figure [5] (a) shows the probability for an inefficient detector to detect two photons 
in the idler mode. The corresponding fidelities and purities have been shown in figure 
[5] (b) and do not vary as a function of \ an d V- 




Figure 5. (a) The probability of detecting two photons in the idler mode for (top to 
bottom): no filter; 07 = 3xl0 12 s _1 ; a f = 2xl0 12 s _1 ; 07 = lxlO 12 s _1 . Note that 
the probability is plotted on a log scale, (b) The purity and fidelity of the signal state 
with an ideal Fock state. 



7. Physical example II - group velocity matching 

In this section, we will examine particular phase matching conditions which result in 
a less entangled JSA, and therefore a more pure heralded Fock state. It is common 
to approximate the phasematching function as <f>(ui,u) s ) = exp(— / ~fL 2 Ak 2 /A) where 
7 w 0.193. By making this approximation, we can write the JSA as 

f{u h u s ) oc exp ^ — j exp ^ J . (87) 
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In order to make equation (87) separable, we require all "cross-terms", i.e. terms which 



contain products of Ui and u s to vanish. This occurs when the condition 
2 

+ lL 2 (k' s -k')(k' i -k'„ 







a- 



(88) 



is met, yielding a JSA of the form f{uji,u s ) oc fi(ui)f s (u s ) [BJ. One way to satisfy the 
condition in equation (88) is to set k' p = {k' s + fcQ/2, which results in the following 
condition for the length of the waveguide, as a function of the pump width: 



L = 1/^87^-^. (89) 
These conditions generate a symmetric JSA, where both signal and idler modes have 



equal widths. Alternatively, rearranging equation (88) as follows 
4 



kt 



+ J<rL(K - k' ] 



0. 



(90) 

k' we can also 



we can see that by making L « a -1 , i.e. L —>■ oo and setting k' p 
obtain a separable JSA. These conditions generate an asymmetric JSA (6j. 

We have made use of the Gaussian approximation for $(ui,u a ) to obtain the 
conditions for separability, however we will now input these conditions into the original 
sine form of the function. This analysis will not result in completely pure states being 
generated, however it should correspond more closely to experimental observations. 



7.1. Symmetric JSA 

In order to meet the extended phase matching conditions for a symmetric, separable 
JSA, we again model a type II PP-KTP waveguide, now of length L = 24.2 mm and 
a periodicity of A = 68.4 /mi, pumped with a 788 nm laser with a 0.7nm FWHM 
(dp = 0.9 xlO 12 s _1 ) which down converts to 1576 nm in the signal and idler modes. 
Unless stated otherwise, the results in this section were obtained using an 800 x 800 
grid, ranging over 0.06 x 10 15 s -1 , centered around u>i = u s = ji. 

Figure [6] shows the JSA and the corresponding Schmidt values and modes. Notice 
in the Schmidt decomposition that the first mode is much more dominant than it was 
in section [6} The entropy of entanglement for this JSA is E = 0.88. If we had used 
the Gaussian approximation for <3>, the JSA would decompose into one pair of Schmidt 
modes and the entropy of entanglement would be E = 0. In such a case, the four-photon 
term would would consist only of two-photon Fock states. 

As an intuitive guide to why the above conditions generate the given JSA, notice 
that varying the parameter Ak has the effect of changing the gradient of the phase 
matching function $(ui,uj s ) (see figure [6] (b)), rotating it around tUi = tu s = /x while 
changing the parameter L, alters the width of the phase matching function. The 
goal is to pick Ak, and therefore k' p , and L such that the phase matching function 
is perpendicular, and of equal width, to the pump function. 

As in the previous section, we present results for the probability, purity, g^ and 
fidelity, for the heralding of one and two photon Fock states. For a realistic JSA, 
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Figure 6. (a) Gaussian pump function a(u>i + u> s ) with a 0.7nm FWHM at fi p = 788 
nm. (b) Phase matching function <fr(uJi, lo s ) for waveguide of length L = 24.2 mm and a 
periodicity of A = 68.4 /xm. (c) The resulting JSA f(u>i, ui s ) = ce(uii+ui s )Q(LL>i, ui s ). The 
JSA has been plotted as a function of the frequency, however corresponding values for 
the wavelength have been included, (d) Schmidt numbers and modes for the JSA (top 
to bottom): the first 20 Schmidt numbers &&; the first 5 Schmidt modes £fc(o> s ) for the 
signal state; the first 5 Schmidt modes Cfc( w i) f° r the ilder state, as well as a Gaussian 
filter function of width cry = lxlO 12 s _1 (dashed line); the filtered Schmidt modes 
T(uJi)(k(oJi) for the idler state; the othogonaliscd idler modes (f>j(uii); the diagonalised 
signal modes T m (u) s ) 



manipulating the phase-matching conditions can result in high purity of the heralded 
state, however, it doesn't reach unity. From figure [6] (c), it can be seen that the outer 
lobes contribute to the spectral correlations and perhaps it is possible to increase the 
purity of the heralded state by filtering them out. Therefore, we will again compare 
results for: an unfiltered idler state; an idler state filtered with a Gaussian filter 
T(uji) = exp(— (fMf — Ui)/2aJ), of various widths Of and centered at the central idler 
frequency; as well as the limiting case where T(uji) = 5(uii — fif)- 

7.1.1. Generating single-photon Fock states Figure [7] (a) shows the probability of 
detecting a single photon in the signal mode. Since most of the photons will have 
spectral distributions within the filter width, we do not see a very big drop in the 
probability, when filtering. 
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Figure 7. (a) The probability of detecting a single photon in the idler mode for (top 
to bottom): no filter; 07 = 3xl0 12 s" 1 ; 07 = 2xl0 12 s" 1 ; 07 = lxlO 12 s" 1 . Note 
that the probability is plotted on a log scale, (b) The fidelity of the signal state with 
an ideal Fock state for (top to bottom at x = an d V = 0): 07 = 0; 07 = lxlO 12 s _1 ; 
07 = 2xl0 12 s _1 ; 07 = 3xl0 12 s" 1 ; no filter, (c) The g (2 ~> of the signal state for (top 
to bottom): 07 = 0; 07 = lxlO 12 s _1 ; cr/ = 2xl0 12 s _1 ; 07 = 3xl0 12 s _1 ; no filter. 
Note the change in axis orientation, (d) The purity of the signal state for (top to 
bottom at \ = and r\ = 0): 07 = 0; 07 — lxlO 12 s _1 ; af = 2xl0 12 s _1 ; no filter. 

Figure [7] (b) shows the fidelity. We distinguish between these surfaces by referring 
to their values at the point \ = and rj = 0. Immediately we can see that the fidelity is 
much higher than in section |6j In the region of interest, filtering the idler mode increases 
the fidelity of the signal state with a single photon. Decreasing the filter width can be 
detrimental to the fidelity when the nonlinearity and the detector efficiency are high. 

Figure[7](c) shows the g^ and figure[7](d) shows the purity of the state in the signal 
mode. Again, we will distinguish between these surfaces by referring to their values at 
the point x — and r] = 0. Results for the filtered case were computed using a 600 x 600 
grid, ranging over 0.06 x 10 15 s _1 , centered around u>i — uj s — /1, and truncating bk with 
values below 10~ 2 . 

Figure [1] (b) shows the probability of success, and required nonlinearity, for a 
number of filter widths, in order to achieve a fidelity of F = 0.95, using a detector 
with efficiency rj = 0.5. Note that while the probability of success is greatly enhanced 
by using a source engineered state, the required pump powers are nearly the same. 
Also notice the "flat" region, where the probability does not change much, between 
07 = 0.6 x 10 12 s _1 and 0.7 x 10 12 s _1 . This corresponds to the "dark" region between 
the lobes on the JSA. Over this region, we do not expect much change in the flux. 
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7.1.2. Generating 2-photon Fock states Figure [8] (a), shows the probability for an 
inefficient detector to detect two photons in the idler mode. The corresponding fidelities 
and purities have been shown in figure [8] (b). The fidelity and purity do not vary as a 
function of x an d t). 




Figure 8. (a) The probability of detecting two photons in the idler mode for (top to 
bottom): no filter; 07 = 3xl0 12 s _1 ; a s = 2xl0 12 s _1 ; 07 = lxlO 12 s -1 . Note that 
the probability is plotted on a log scale, (b) The purity and fidelity of the signal state 
with an ideal Fock state. 



7.2. Asymmetric JSA 

In order to meet the extended phase matching conditions for an asymmetric separable 
JSA, we again analyse a type II PP-KTP waveguide, of length L = 80 mm and 
a periodicity of A = 232 jum, pumped with a 1.93 jum laser with a 3nm FWHM 
(cr p = 0.64xl0 12 s -1 ) which down converts to 3.85 /im in the signal and idler modes. 
We note that single-photon detection is not particularly practical at this wavelength, 
however for consistency, we have chosen to use a PP-KTP waveguide throughout 
this paper. The same JSA can be achieved in different systems, at more practical 
wavelengths. See, for example, Mosley et al. [14]. Unless stated otherwise, the results 
in this section were obtained using an 800 x 800 grid, ranging over 8xl0 12 s _1 , centered 
around u)i = u s = fi. 

Figure [9] shows the JSA and the corresponding Schmidt values and modes. Notice 
in the Schmidt decomposition that the first mode is even more dominant than in the 
symmetric case. The entropy of entanglement for this JSA is E = 0.37. 

Setting k' p = k' s generates a vertical phase matching function. As long as the 
waveguide is sufficiently long, and therefore, the width of the phase matching function 
sufficiently thin, and the pump is sufficiently wide, the result will be a vertical, almost 
elliptical and very thin JSA. 

As in the previous section, we present results for the probability, purity, g^ and 
fidelity, for the heralding of one and two photon Fock states. Although it is possible 
to achieve purities arbitrarily close to unity by increasing the length of the waveguide, 
the vertical orientation of the JSA places it in a unique position to take advantage of 
spectral filtering. We will once again compare results for: an unfiltered idler state; 
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Figure 9. (a) Gaussian pump function a(uji + uj s ) with a 3nm FWHM at fi p = 
1.93 fim. (b) Phase matching function Q(oji,uj s ) for waveguide of length L = 80 mm 
and a periodicity of A = 232 /im. (c) The resulting JSA f(u>i, uj s ) — a(u)i+uj s )Q(u>i, u> s ). 
The JSA has been plotted as a function of the frequency, however corresponding values 
for the wavelength have been included, (d) Schmidt numbers and modes for the JSA 
(top to bottom): the first 20 Schmidt numbers bk', the first 5 Schmidt modes £fc(w s ) 
for the signal state; the first 5 Schmidt modes CfcC^i) f° r the ilder state, as well as 
a Gaussian filter function of width <jf — 100 xlO 9 s _1 (dashed line); the filtered 
Schmidt modes T(yii)C,k(^i) for the idler state; the othogonalised idler modes <j)j{u>i); 
the diagonalised signal modes T m {ui s ) 



an idler state filtered with a Gaussian filter T(ui) = exp(— — Ui)/2a^), of various 
widths o~f and centered at the central idler frequency; as well as the limiting case where 
T(ui) = 5(u)i - /if). 



7.2.1. Generating single-photon Fock states Figure 10 (a) shows the probability of 



detecting a single photon in the signal mode. Since most of the photons will have 
spectral distributions within the filter width, we do not see a very big drop in the 
probability with filtering, until the filter is so narrow that it cuts into the central lobe. 



Figure 10 (b) shows the fidelity. We will distinguish between these surfaces by 
referring to their values at % — 0.5 and 77 = 1. Again, filtering the idler mode increases 
the fidelity of the signal state with a single photon in the low-x and low-r? regimes. 

Figure 10 (c) shows the g^ and figure 10 (d) shows the purity of the state in the 
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Figure 10. (a) The probability of detecting a single photon in the idler mode for (top 
to bottom): no filter; 07 = 150xl0 9 s" 1 ; a s = lOOxlO 9 s _1 ; 07 = 50xl0 9 s^ 1 . Note 
that the probability is plotted on a log scale, (b) The fidelity of the signal state with an 
ideal Fock state for (top to bottom at X — 0-5 and r\ = 1): no filter; 07 = 150xl0 9 s _1 ; 
07 = 100 xlO 9 s" 1 ; 07 = 50xl0 9 s" 1 ; 07 = 0. (c) The g^ of the signal state for (top 
to bottom): 07 = 0; 07 = 50xl0 9 s" 1 ; 07 = 100x10 s s" 1 ; 07 = 150xl0 9 s _1 ; no 
filter. Note the change in axis orientation, (d) The purity of the signal state for (top 
to bottom at x = 0-5 and 77 = 1): no filter; 07 = 150xl0 9 s _1 ; 07 = 100xl0 9 s -1 ; 
a f = 50xl0 9 s" 1 ; 07 = 0. 

signal mode. Again, we will distinguish between these surfaces by referring to their 
values at the point \ — 0-5 and 77 = 1. Results for the filtered case were computed using 
a 600 x 600 grid, ranging over 8xl0 12 s -1 , centered around Ui = uj s = /j, and truncating 
bk with values below 10~ 2 . 

Figure|4](c) shows the probability of success, and required nonlinearity, for a number 
of filter widths, in order to achieve a fidelity of F — 0.95, using a heralding detector 
with efficiency 77 = 0.5. Note that, as with the symmetric case, while the probability 
of success is greatly enhanced by using a source engineered state, the required pump 
powers are nearly the same. Also notice the "flat" region, where the probability does 
not change much, between 07 = 0.3 x 10 12 s _1 and 0.35 x 10 12 s _1 . This corresponds 
to the "dark" region between the lobes on the JSA. Over this region, we do not expect 
much change in the flux. 
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7.3. Generating 2-photon Fock states 



Figure 11 (a), represents the probability of detecting two photons in the idler mode. 
The corresponding fidelities and purities have been shown in figure [8] (b). They do not 
vary as a function of x an d fj. 
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Figure 11. (a) The probability of detecting two photons in the idler mode for (top 
to bottom): no filter; 07 = 150xl0 9 s -1 ; 07 = lOOxlO 9 s" 1 ; 07 = 50xl0 9 s -1 . Note 
that the probability is plotted on a log scale, (b) The purity and fidelity of the signal 
state with an ideal Fock state. 



8. Discussion 

We have calculated the spectrally entangled output state of a parametric down converter 
to second order in photon number, with the goal of generating heralded one- and two- 
photon Fock states in one spatial mode (signal), conditional on the detection of one or 
two photons in the other spatial mode (idler). We have presented analytical expressions 
for the heralded state after the idler mode is spectrally filtered using a Gaussian filter and 
detected with an inefficient detector. The heralded signal state was then characterised 
by its g^ and purity. In addition, we calculated the fidelity of the heralded state with 
the desired ideal Fock state. 

As a physical example, we modeled a type II PP-KTP waveguide, pumped by lasers 
at wavelengths of 400 nm, 788 nm and 1.93 fim. We found that in the first example, 
where no effort was made to perform any extended phase matching conditions, the 
results were states with very low purity. After strong spectral filtering, Fock states with 
arbitrarily high purity could be achieved, however at very low probabilities of success. 
To achieve a fidelity of F = 0.95 for a single-photon state, using a heralding detector 
with efficiency 77 = 0.5, the probability of success would be on the order of 10 -4 . 

The latter two examples, where extended phase matching conditions were fulfilled, 
resulted in much higher purity states, however, some additional filtering was still 
required to achieve very high purity states. Both the symmetric and asymmetric 
examples were able to achieve a fidelity, with a single-photon Fock state, of F = 0.95, 
using a heralding detector with efficiency r\ = 0.5, with probabilities of success on the 
order of 10~ 2 . High purity two-photon Fock states were also possible. While results were 
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comparable for the symmetric and asymmetric examples, this was due to our choice of 
physical parameters. The asymmetric case would be able to achieve higher fidelities, 
with no filtering, by choosing a longer waveguide. 
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